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Abstract 

We discuss linear in momenta Poisson structure for the generalized nonholonomic Chap- 
lygin sphere problem and prove that it is non-trivial deformation of the canonical Poisson 
structure on e* (3) . 

1 Introduction 

Let us consider a rolling of dynamically asymmetric and balanced spherical rigid body, the 
so-called Chaplygin ball, over an absolutely rough fixed sphere with radius a pQ . At a — > oo one 
gets a Chaplygin problem on a non-homogeneous sphere rolling over a horizontal plane without 
slipping g]. 

Since slipping at the contact points is absent, its velocity vanishes and we have the following 
nonholonomic constraint 

v + uj x r = . (1) 

Here uj and v are the angular velocity and velocity of the center of mass of the ball, r is the 
vector joining the center of mass with the contact point and x means the vector product in R 3 . 
Mass, inertia tensor and radius of the rolling ball will be denoted by to, I — diag(Ji, I2, 13) and 
b, respectively. 

According to [I] , the angular momentum M of the ball with respect to the contact point 
with the sphere is equal to 

M = (I + dE) uj - d(7, w)7 , d = mb 2 . (2) 

Here 7 is the unit normal vector to the fixed sphere at the contact point, E is the unit matrix 
and (., .) means the standard scalar product in M 3 . All the vectors are expressed in the so-called 
body frame, which is firmly attached to the ball, its origin is located at the center of mass of 
the body, and its axes coincide with the principal inertia axes of the body. 

After elimination of the Lagrangian multiplier according to [T], one gets the following 
reduced equations of motion 

M = M x uj , 7 = K7 x uj , where K = . (3) 

a + b 

These equations possess three integrals of motion 

Hi = (M,«), H 2 = (M,M), d = (7,7), (4) 

and invariant measure 



where 



/1 = g- 1 ( 7 ) d 7 dAf , g( 7 ) = Vl-d(7,A 7 ), (5) 
(I + dE)- X =| ^ 
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At k — il one more integral of motion exists 

h 

C 2 = ( 7 ,BM), B=| b 2 I =trA- 1 + ( K -l)A- i . (6) 



6. 
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It is easy to see, that in the Chaplygin case a — > oo we have re = 1. 

At re = — 1 we have generalised Chaplygin sphere problem or so-called Borisov-Mamaev- 
Fedorov system, see PQ and [3]. 



2 The Poisson brackets 

At re = ±1 six equations of motion ([3]) possess four integrals of motion and an invariant measure. 
Thus, by the Euler-Jacobi theorem, they are integrable in quadratures. It allows us to suppose 
that common level surfaces of integrals form a direct sum of symplectic and lagrangian foliations 
of dual dynamical system which is hamiltonian with respect to the Poisson bivector P, so that 

[P,P]=0, PdCi, 2 = 0, (PdH 1 ,dH 2 ) = {H 1 ,H 2 } = 0. (1) 

Here [., .] is the Schouten bracket. In fact, we suppose that the Euler-Jacobi integrability of 
non-Hamiltonian system (|3j is equivalent to the Liouville integrability of the dual Hamiltonian 
dynamical system with the same integrals of motion, see [5] . 

The first equation in (TIJ guaranties that P is a Poisson bivector. In the second equation we 
define two Casimir elements C\,i of P and assume that rankP = 4. It is a necessary condition 
because by fixing its values one gets four dimensional symplectic phase space of the desired 
Hamiltonian system. The third equation provides that integrals Hi 2 are in involution with 
respect to the Poisson bracket associated with P and, therefore, that they form a lagrangian 
foliation. 

In order to compare Poisson structures at re = ±1 we briefly remind some known facts 
about linear in momenta M solutions P of (fTJ associated with the Chaplygin sphere problem 
at re = 1 following to [3 El El [7]. 



2.1 Chaplygin sphere, k = 1 . 

According to [2] , integrals of motion (gJE]) are in involution with respect to the following Poisson 
brackets 

/ d(M At) \ 

{Mi, M } g = e ijk gM k - '- " lk , {M i)7i } 9 = e ijk gjk , lj} a = 0, (2) 



where e^k is a totally skew-symmetric tensor. These brackets have the necessary Casimir 
functions C\^- 

In variables x = (71, 72, 73, M\, M%, M3) initial equations of motion (|3]) have the form 

^=X h =g- l {H,x k } g , where H = ^ . (3) 
After a change of time 

At -> gdt (4) 

these equations becomes Hamiltonian equations with respect to the Poisson brackets ([2]). It 
means that initial non-Hamiltonian vector field X is the conformally Hamiltonian vector field 

X = g- 1 (x) X , where X = P 9 dH . 

The Poisson brackets ([2]) can be easily obtained via trivial deformations of the canonical 
Poisson brackets and the standard momentum map theory. Namely, let Q be a n-dimensional 
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smooth manifold. Its cotangent bundle T*Q is naturally endowed with the Liouville 1-form 9 
and symplectic 2-form 57 = d9, whose associated Poisson bivector will be denoted with P. In 
local symplectic coordinates on T*Q 

z = (q,p) = (qi,---,qn,pi,---,p n ) 

they have the following local expressions 

6 = pi&qi + ... Pndq n , fi = d# = pi A Qi H p n A q n . 

Let us substitute the scaling momenta 

Pk^g(q)pk k = l,...,n, (5) 
into the Liouville and symplectic forms 

9 g = g(q) (pidqi + ■ ■ -p^Aq^j , ft -> n g = de g . 

The corresponding Poisson bracket 

tii , Qj } g = , {qt , pj } g = gS.j , { Pl , Pj } g = d j gp l - d l gp j , (6) 

is a trivial deformation of canonical Poisson bracket 

{q h qj} = 0, {quPj} = hj , {Pi,Pj} = 

in the Poisson-Lichnerowicz cohomology |S]. Here dk = d/dqk- 

Now let us identify Q with a two dimensional sphere S 2 embedded into K 3 , so that q^ — ji, 
i = 1, 2, 3. The standard momentum map 

4> : (p, 7) eTS 2 4 (M, 7) G e*(3) = so(3) k M 3 , 

defined by the vector product 

M = jxp, (7) 
maps our trivial deformation (J6j) into the following Poisson brackets on the Lie algebra e*(3) 
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{Mi , Mj } = s ijk g(-y)M k + lk ^ M m d m g(^) , (8) 



rn— 1 



{ A ^ » 7? } ff = £ yft 3(7) 7fc , {t* , 7j } g = • 
As above, here we use an abbreviation d m — d/dj m . 

Proposition 1 If we identify 5(7) with g("f) then the Poisson brackets £5P coincide with 
the Poisson brackets (0). 

So, the Poisson brackets © are trivial deformations of canonical ones. Consequently, according 
to [7J, change of variables 

L 1 =^(m 1 -P^(i + ^ ' C71 



(7, 7) V v ) ) 7l 2 + 75 

2 s 



3 



where 




b=(j,M), c=(7,i) and v = j 2 + j 2 - d(j, 7)(ai7i + 0272) > 

allows us to reduce the deformed Poisson brackets © to the canonical Lie-Poisson brackets on 
the Lie algebra e*(3) 

{L l ,Lj} = Ei-jkLk , {L t ,Jj} = Stjkjk , {li , 7j } = . (10) 

So, we can prove that the original nonholonomic Chaplygin system is trajectory equivalent to 
the dual integrable dynamical system on two-dimensional sphere S 2 , which is a Hamiltonian 
system with respect to the canonical Lie-Poisson brackets (fTU|) . see details in [7J. 

2.2 Generalized Chaplygin sphere, k = —1 . 

Now let us compare known Poisson structure at k — 1 with the new Poisson structure obtained 
for the case k = — 1. It is easy to see that at k = 1 bivector P g associated with the Poisson 
brackets @ has the form 

^ 9 = g(_r T m)-^ 1(m ' A7) (o r)' 

where 

73 -72 \ 
-73 71 , M 
72 -71 / 

Of course, canonical Poisson bivector on e*(3) 

P ={-T- m)' ( 12 ) 
is compatible with its trivial deformation P g so that 

[P,Pg]=0. 

The Poisson bivector P& for the generalized Chaplygin ball rolling over the sphere, albeit on 
the similar form, has completely another properties. 

Proposition 2 At k = — 1 integrals of motion (0|6p are in involution with respect to the 
Poisson brackets defined by the following Poisson bivector 

n = g I I +g- 1 (2rf(7,7)-trB)^ ~ ) , (13) 

which is just one linear in M solution of the eguations fl]). Matrix T depends only on 7 

(7,7)E-C-g 

where E is a unit matrix, 

7i 7i72 7i 73 
7271 72 7273 I and 
737i 7372 73 
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Entries of the other two matrices are equal to 

blMk 
2d 



My = -eijk a k -f k - (7,7)6 fc M fc + 



f « - ' + +ff)fr + H) ( {k + bj) ak + {bk ~ bi){hk ~ W7fe ) 

here 

a k = (c 2 + 6 fc (7, M)) , C 2 = (7, BM) . 
The corresponding Poisson brackets look like 
{H,lj}b = 0, 

M7l + 7l) + fo 272 - 

^2(7? +7l) +&37I - ~2^- 
{M 2 , 7 2 }b = g(6 3 ~ ^1)717273 , {M 2 , 7i} fc = -g73 (b%7i + hill + 7a) ~ 

{Af 2 , 73 } fc = g 7 i ( 61(7? + 7I) + 637I 



2d 

i M:;. ~.:\}h = ~ '>2 I" l7j7; - ! -U;. " I i /. = JO 2 ( &l7l + ^2 (tI + 73 ) ~ 



{M 3) 72 } fc = -g 7l (^ (7l 2 + 7 f) + 6 27 | " ^ ) 

and 

{M 1 ,M 2 } b = g(-a373 + (7,7)&3M 3 -%^) - g~ 1 7 3 (2d(7, 7) - trB) x 



0363 1 (b 3 -b 2 )(b 3 ~b 1 )b 3 -y 3 M 3 \ 

(b 3 +b 2 )(b 3 +b 1 ) (b 3 +b 2 )(b 1 +b 3 )(b 2 +b 1 ) J 



Brackets {Mi,M$}b and {M 2 , Ms}b have the same form as {M\, M 2 }b and, therefore, we omit 
their explicit expressions. 

If C 2 — there are many other linear in momenta M solutions of the equations ([TJ 
associated with known variables of separation, see details in [BJ. 

At C 2 7^ using linear in momenta Poisson brackets {-,-}b we can rewrite equations of 
motion ([3]) in the following form 

^ = X k = g^{H u x k } b + g^{H 2 ,x k } b , (14) 

where 
and 

s( 7 ) =4d 2 (7,7)(7,B 7 )-2d((EtrB-B)7,B 7 ) + detB . (15) 

It is easy to see that at k = — 1 equations of motion have a more complicated from in comparison 
with the original Chaplygin problem Q at k = 1. 
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Proposition 3 At k = — 1 the initial non-Hamiltonian vector field X hS[14\j is a sum of two 
conformally Hamiltonian vector fields 

x= g r 1 ii+g 2 - 1 x 2 , 

where X\ }2 are hamiltonian vector fields associated with two commuting integrals of motion 

Xx = P b dffi , X 2 = P b dH 2 , {H u H 2 } b = . 
On the other hand, this non-Hamiltonian vector field is conformally Hamiltonian vector field 

X = gs 1 X 3 , where X 3 = P b dH 3 and g 3 ( 7 ) = lilMl) , ( 16 ) 
but with respect to another Hamiltonian 

H 3 = (2d( 7 , 7) " trB)#i + 2H 2 , (17) 
which is an integral of motion of f3J) without any distinguished physical meaning. 

At K = 1 the similar linear combination H 3 = H 2 — dH\ coincides with the Hamiltonian of 
the Veselova system, which is equivalent to the original Chaplygin ball at the special choice of 
parameters a% [7]. 

Now let us discuss the main difference between bivectors P g and Pb- 

Proposition 4 Bivector P b \1S\) is nontrivial deformation of the canonical Poisson bivector P 
ftlty) on the Lie algebra e*(3). 

In contrast with bivector P g (TIT]) bivector P b (fT3"|) is incompatible with the canonical Poisson 
bivector P (JT2J) on e*(3) because 

[P, P g ] = , whereas [P, P b ] ± . 

As sequence bivector P b can not be trivial deformation of P. 

Nevertheless, bivector P b is nontrivial deformation of canonical bivector P, because there 
is change of variables 

L i = nTlHTTT 2^ f«i (6172^1 ~ b 2ll M 2 ) + ftMa + + -ff-^ , 

g( 7 )s(7)6i6 2 (7i + 7 2 ) V 7i + 1 2 ) 7i + T 2 

L2 = < \ ru 1. ? 2 t ^ f « 2 (&i7 2 Mi - 6271M2) + feM 3 + 57 2 7 ? fe y ) + -riS ■ 

g(7)a(7)&i fe(7i + 7a) V 7i + 72 / 7i + 7 2 

is = , w vt 7 , 2 , 2 , (a 3 (fei72Mi - 6 27 iAf 2 ) + 3 M 3 - 6/1(7)) , 
g(7)s(7)6i6 2 (7i + 7 2 ) v / 

b = (S 7 ,M), c=( 7 ,i), 
which allows us to reduces this bivector to canonical one. Here 5(7) is given by (|15p . 

ai = 2dj 2 (2dbijf + 2db 2 (j 2 + 7 2 ) - bib 2 ) , a 3 = 4d 2 7i 7^3(61 - 6 2 ) , 

q 2 = -2rf 7l (2d6i( 7l 2 + 7 2 ) + 2d6 27 2 - 6162) , 

and 

ft = -2d7 1 7 3 (2rf(7 1 2 (6ife 2 - 6i6 3 + 6 2 6 3 ) + 62(617? + 6373)) - 616263) , 

(3 2 = -2d7 2 7 3 (2d(7 2 (6i63 + b x b 2 - b 2 b 3 ) + h(b 2 ^ + 6373)) - 6i6 2 6 3 ) , 

03 = 4d 2 6 3 73 2 (6 2 7i 2 + 617I) + 2d6 1 6 2 ( 7l 2 + 7 2 2 ) (2d(7? + 7 2 2 ) - 63) . 
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For the brevity we omit the explicit expressions for the function ^.(7), which is a solution of the 
differential equations {i.; , Lj } fc = EijkLk ■ 

Applying this transformation at c — (7, L) = to the Hamilton function (|17p one gets 
integrable dynamical system on the two-dimensional sphere, which is a standard Hamiltonian 
system with respect to canonical Poisson brackets on T*S 2 . The generalised nonholonomic 
Chaplygin sphere is trajectory equivalent to this Hamiltonian system up to change of time 
defined by (fTffjl . 
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